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Abstract 



t3 

We propose an atom-superposition-based method for simulating spin-polarized scanning tun- 

o : 

O , neling microscopy (SP-STM) from first principles. Our approach provides bias dependent STM 



images in high spatial resolution, with the capability of using either constant current or constant 
height modes of STM. In addition, topographic and magnetic contributions can clearly be dis- 



00 . 

I/-) , tinguished, which are directly comparable to results of SP-STM experiments in the differential 



magnetic mode. Advantages of the proposed method are that it is computationally cheap, it is 
easy to parallelize, and it can employ the results of any ab initio electronic structure code. Its 
capabilities are illustrated for the prototype frustrated hexagonal antiferromagnetic system, Cr 
k> \ monolayer on Ag(lll) in a noncollinear magnetic 120° Neel state. We show evidence that the 

magnetic contrast is sensitive to the tip electronic structure, and this contrast can be reversed 
depending on the bias voltage. 
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I. INTRODUCTION 
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Research on magnetic systems has intensified in the last decade with the u 
to produce magnetic data storage devices with ultrahigh information density 
can be achieved by reducing the size of the information storage units going down to the 
nanoscale or even to single atoms 3]. Detecting and manipulating spins [4J with high 
accuracy on the atomic scale is, thus, essential for future technological applications. Spin- 
polarized scanning tunneling microscopy (SP-STM) |5| is one of the main tools for studying 
magnetism at the atomic scale. Recent experimental advances using this technique allow the 
investigation of complex magnetic structures (frustrated antiferromagnets, spin spirals, etc.) 
6|— llljl . Considering such structures in reduced dimensions, their magnetic ground state can 
be determined [7|, |ll| and the nature of magnetic interactions can be studied by theoretical 



means, e.g. from first principles (12(, or applying a multiscale approach [131 ] . However, a 
proper validation of the proposed ground state spin structures demands a method, which 
is capable to directly compare them to experimental observations. This can be done by 
SP-STM simulations. 

Our motivation was to construct a computationally cheap and user-friendly, yet reliable 
model for simulating SP-STM. Here, we propose an efficient method based on the spin- 



polarized Tersoff-Hamann model |l4j and the atom superposition approach 15l-ll7|. Our 
model goes beyond the work of Heinze (l5|, and considers in particular, 

(1) bias voltage, for simulating bias dependent physical properties, 

(2) energy dependence of the vacuum decay of electron states, and 

(3) energy dependence of atomic local spin quantization axes. 

The main advance of our tunneling model is the inclusion of the tip electronic structure, 
which is neglected in Refs. |14l.ll5|. and it enables to study tip effects on the SP-STM images. 
First, we determine energy dependent virtual differential conductance (dl/dU) quantities 
on a three-dimensional fine real space grid from electron local density of states (LDOS). 
Integrating the differential current (dl) contributions in an energy window in accordance 
with the applied bias voltage, we obtain a three-dimensional current map on the same 
grid, from which constant current or constant height images can be extracted. Calculating 
differential current first and then tunneling current by integration over energies proved to 
be numerically more stable than the opposite (numerical differentiation of tunneling current 



with respect to energy) 18[. Furthermore, bias dependent apparent barrier height can be 
included in our approach. 

The paper is organized as follows: The theoretical model of SP-STM ispresented in 
section [Til where we explicitly point out extensions to the Heinze model [l5[ and other 



atom-superposition-based approaches 16|, |l7j. As an application, we investigate the frus- 



trated hexagonal antiferromagnet, one monolayer (ML) Cr on Ag(lll) in section [TTTl We 
study two magnetic chiralities for its noncollinear Neel state and determine the energetically 
favored magnetic structure. Moreover, we compare electronic structures of Cr obtained 
from collinear and noncollinear calculations. By simulating SP-STM images we are able to 
investigate magnetic contrast depending on tip electronic structure and bias voltage. Our 
conclusions are found in section IIVI 

II. THEORETICAL MODEL OF SP-STM 



In SP-STM the total tunnelin g cu rrent can 



Io, and a spin-polarized part, Ip 



14 



15 
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3e written as a sum of a non-spin-polarized, 



h + h 



(1) 



This formula is generally valid for either collinear or noncollinear surface and tip spin struc- 
tures. Iq and Ip can be calculated at different levels of approximation for the tunneling 
current. It could, in principle, be implemented within the multiple scattering framework 
20[. While the perturbation approach has been used by Hofer and Fisher 19] for collinear 
surface and tip spin structures with an arbitrary angle between their spin quantization axes, 



the most commonly used method is based on the Tersoff-Hamann model 21, 221. Wortmann 

n 

et al. [14J introduced its spin-polarized version applicable to complex noncollinear surface 
spin structures. Heinze 15[ combined this with the atom superposition method jig, Il7| . 
Note that in the following we denote I and Ip by Itopo (topographic current) and Imagn 
(magnetic current), respectively. 



Here, based on the work of Heinze 15J, we propose a hybrid model, which uses essentially 
the Tersoff-Hamann formalism but we do not restrict the tip electron density of states (DOS) 
to be constant in energy. This means that different tip models 23[ and their effect on 
tunneling properties can be investigated. The only requirement for our present formalism is 



that we assume that electrons tunnel through one tip apex atom. Since the tip electronic 
structure is explicitly included in our method via the projected DOS onto this apex atom, 



Eq.(2) of Ref. [15| needs to be reconsidered. Our strategy is to determine differential currents 
first and then perform an energy integral in a window according to the applied bias voltage 
(V) in order to arrive at the tunneling current. 

Let us define the following position- and energy- dependent density matrices in spin space 
for tip (T) and sample (S), respectively, 

EpgfaE) = n Ti5 (r,E)I + m Ti5 (r,E)^ (2) 

^t,s(^; E) + m T s (r, E) m T s (f, E) — im TS (r, E) 
m T s (r , E) + im T s (f,E) n TtS (r , E) — m T s (f,E) 

Here, J is the 2x2 unit matrix, W is the Pauli matrix vector, while ut{Rtip,E) and 
rnxiRriPjE) are the charge and magnetization DOS projected to the tip apex atom. On 
the other hand, ns{f, E) and ms(r, E) are the charge and magnetization DOS of the sample 
surface at position f. They can be obtained from the corresponding density matrix as 

n T ,s(r,E) = ±Tr(p s (r,E)y (3) 

m T)S (r, E) = l -Tr (p^f, E)EJ , (4) 

where the Trace is performed in the spin space. At the tip position Rtip above the surface 
we obtain the charge and magnetization electron local density of states (LDOS) of the 
surface in vacuum, us(Rtip,E) and ms(RTip,E), respectively. Combining the sample and 
tip density matrices at Rtip, a modified LDOS can be defined as 

LDOS(R TIPl E) = AE^Tr (^(Rtip, E)^(R tip , £)) (5) 

= AE (n s (R T iP, E)n T (R TIP , E) + fn s (R TIP , E)m T (R TIP , E)) , 

which, in fact, combines the vacuum LDOS of the surface and the projected DOS of th e tip 



apex atom. This formula is consistent with the spin-polarized Tersoff-Hamann model 14|, 



except the fact that it explicitly includes the electronic structure of the tip apex. Here, AE 



ensures that the LDOS is correctly measured in units of (eV)~ l . Note that in Ref. 15| the 
tunneling current was proportional to a dimensionless LDOS at the Fermi level. 

The vacuum LDOS of the surface can be approximated by a superposition of decaying 
atomic electron states. Following this, we consider the position dependence of the sample 



density matrices as p (R a , E) with R a the position vector of the ath sample surface atom, in 
order to allow different chemical or magnetic properties for these atoms. n^{E) = ns{R a , E) 
and rfTg(E) = fns(R a ,E) now denote charge and magnetization DOS projected to the ath 
surface atom, respectively. It has to be noted that chemical differences between surface 



atoms were not taken into account in Ref . 15J . A tunneling transition between the tip apex 



and the ath surface atom at energy E can be represented as the Trace of the multiplied 
density matrices, similarly to Eq.([5]). This is the energetic ingredient for the tunneling 
transition. Apart from this, the transmission coefficient through a potential barrier between 
the ath surface atom and the tip apex has to be included in the tunneling model, in general 
we denote it by T(E, V, d a ). It has energy and bias dependence, and contains geometry 
information of the three-dimensional tunnel junction via the distance between the tip apex 
and the ath surface atom, 

d a (x,y,z) = \R TIP (x,y,z) -R a \- (6) 

Thus, the modified LDOS at the tip apex position R TIP (x,y, z) and at energy E can 
be approximated as the superposition of individual atomic contributions from the sample 
surface as 

LDOS(x, y, z, E, V) = AE ^ T(E, V, d a (x, y, z))^Tr (^(Ra, E)pJR TIP , £)) . (7) 

a 

The main advantage of using the density matrix formalism is that electronic and spin struc- 
tures calculated either nonrelativistically or relativistically can be treated within the same 
theoretical framework. 

We calculate the above LDOS values at (x,y, z) grid points of a three-dimensional fine 
grid in a finite box above the surface. The image resolution is determined by the density of 
(x, y) grid points. The motivation for using the atomic superposition approximation is, on 
one hand, computational efficiency, since calculating and storing the projected DOS onto 
surface atoms in the magnetic unit cell is computationally much cheaper compared to the 
vacuum LDOS of the surface on a great number of grid points. On the other hand, such 
atom-projected DOS functions are routinely obtained in all ab initio electronic structure 
codes, whereas vacuum LDOS is not always routinely accessible for the average user. 

According to above, the LDOS can be decomposed, similarly to Eq.flT]), as 

LDOS(x, y, z, E, V) = LDOS TOPO (x } y, z, E, V) + LDOS MAGN (x, y, z, E, V), (8) 



and assuming an exponential decay of the electron wavefunctions, the TOPO and MAGN 
terms can be written as 

LDOS TO po(x, y, z, E, V) = AE ^ e- 2 « E > v ^^nr(E)n a s (E) (9) 

Q 

LDOS MAGN (x,y,z,E,V) = AEj2e~ ME ' V)da{x ' v ' z) ™T(E)m a s (E)cos<p a (E). (10) 

a 

Here the sum over a has to be carried out, in principle, over all the surface atoms. Con- 
vergence tests, however, showed that including a relatively small number of atoms in the 



sum provides converged LDOS values [24|. Each surface atom is characterized by a local 
spin quantization axis, eg(E), which can be defined from the sample magnetization DOS 
vector as e$(E) = fffg(E) /\rffg(E)\. In the most general case, these local axes can be energy 
dependent, see Eqs. (TT9l and (120|) . and Table fl] This also holds for the spin quantization axis 
of the tip apex, ct{E). The exponential factor is the transmission probability for electrons 
tunneling between states of atom a on the surface and the tip apex, where k is the vacuum 

nnn 

decay k is treated within the independent-orbital approximation [15l I21L |22[ , which means 
that the same (spherical) decay is used for all type of orbitals, but its energy dependence is 
explicitly considered essentially in the same fashion as in Ref. [25[. Extension of our model 
to take into account orbital dependent vacuum decay following Chen's work [26] is planned 
in the future, which is relevant for a more advanced description of tunneling from/to direc- 
tional orbitals. In the present paper we consider two different ways of calculating k. One is 
inspired by the Tersoff-Hamann model, taking only surface properties into account, 



K(E) = ^2m( ( Ps-\E-E s F \), (11) 

where the electron's mass is m, K is the reduced Planck constant, while (fis and Ep are 
the average electron workfunction and the Fermi energy of the sample surface, respectively. 
The absolute value ensures that the transmission probability is symmetric in the positive 



and negative bias range with the minimum at zero bias [27|, [28|] . This way we circumvent 
the problem of the bias-asymmetric contribution from the sample LDOS to the differential 
conductance [29j. We use this energy dependent vacuum decay for an ideal, electronically 
featureless and maximally spin-polarized tip model. Note that this formula does not have an 
explicit bias dependence. Taking into account the tip apex electronic structure obtained from 
first principles, the more general expression for k is based on the one-dimensional Wentzel- 
Kramers-Brillouin (WKB) approximation assuming an effective rectangular potential barrier 
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between the tip and the surface, 



K (E,V) = ^2m ^ + ^ + eV -\E-E^, (12) 

with <pT being the local electron workfunction of the tip apex, e the elementary charge and 
V the applied bias voltage. This vacuum decay formula is considered for our magnetic Ni tip 
model. The quantity (05 + 0y + eV A )/2 — \E — E§\ is the energy- and bias dependent apparent 
barrier height for tunneling electrons, 4>n(E, V). Empirical or model nonlinear variations of 



a {E, V) with respect to bias voltage 27|, |28( can also be included in our approach. Note 



that in the case of (pT + ^V = <ps, the first expression of n, Eg. (II II) is recovered. The average 
workfunction of the sample surface is calculated from the local electrostatic potential on a 
three-dimensional fine grid, $(x, y, z), as 



max 



{aj^E*<™*>}-^ (13) 



with N x and N y the corresponding number of grid points, and the local workfunction of the 
tip apex is obtained as 

(p T = max{®(x ,y ,z)} - E%, (14) 

z 

with x and y§ lateral coordinates of the tip apex atom, and Ej, the Fermi energy of the tip 
material. 

In the LDOS formula Eq. ffHJj) . ip a (E) is the angle between the spin quantization axes of 
the tip apex and the ath surface atom at energy E. Previously, only the case of energy inde- 
pendent if a has been considered [151 ]. which corresponds to the angle between the directions 
of local magnetic moments of surface atoms and the tip magnetic moment. However, there 
are more possibilities to combine electronic structure data of sample and tip, which may re- 
sult in an energy dependent (p a (E), see Table [D All listed combinations can be investigated 
within our formalism. The combinations framed by black solid lines are considered in the 
present work, while the one with the gray dashed line corresponds to the studied system in 



Ref. 24|. In Eq. ffTOl) . nir r {E) and rrig(E) denote electron magnetization DOS projected to 



the tip apex and the orth surface atom, respectively, in the collinear case, 

m TjS {E) = n\ s {E) - n l TS (E), (15) 



■f and 4- relative to their local spin quantization axes. Similarly, in Eq.(j9]), ht(E) and Tig(E) 
are electron charge DOS projected to the tip apex and the ath surface atom, respectively, 

n T)S {E)=nl tS {E)+n^ s {E). (16) 

The spin- resolved atom-projected DOS (PDOS) quantities, rtj * S (E), are obtained from first 
principles collinear magnetic calculations. For this task any available ab initio electronic 
structure code can be used. This flexibility of the present SP-STM approach is expected 
to be highly advantageous. Spin-resolved PDOS is considered by assuming a Gaussian 
broadening of the peaks at the k-resolved spin-dependent electron energy (Kohn-Sham) 
eigenvalues, £?r'g(k), obtained at zero temperature, as 

4i(*0 = E E ^e-(-4!^)) V I d ^§^§\r), (17) 

* •* atomic volume 

with ipTS v) ^ e spin-dependent electron wavefunctions corresponding to e ] j,g(k) for tip 
(T) and surface (S), respectively, and j the energy band index. The integral over the atomic 
volumes can be performed either in the atomic sphere or within the Bader volume 30j. In 
the present study we use integral over atomic spheres. The Gaussian parameter G could, in 
general, be temperature dependent. In our calculations, we fixed it to a relatively high value 
of 0.1 eV in order to provide smooth n^ S (E) functions. Concerning smoothness of PDOS, 
a high G value counteracts the effect of eventually underrepresented bulk states due to a 
slab geometry, and it is useful if the number of k-points in the Brillouin zone is restricted 
due to computational reasons. 

As Heinze pointed out 15[ , in case of having chemically equivalent surface atoms, the spin 
structure plays a much more dominant role compared to the detailed electronic structure in 
determining the main features of an SP-STM image. This means that SP-STM simulation of 
a known noncollinear spin structure can reasonably be approximated based on the collinear 
electronic structure. We would like to check this statement, therefore, in this paragraph we 
show how to incorporate the fully noncollinear electronic structure into our model. In this 
case the atom-projected charge DOS at energy E is obtained in the following way, 

nr,s{E) = E E ^e-(^-®) 2 / G2 J dV*$(r)¥# (r), (18) 

™ ' atomic volume 

where e^ g (fc) is the set of electron energy (Kohn-Sham) eigenvalues at zero temperature, and 
$ J rs (f) the corresponding spinor electron wavefunctions. The atom-projected magnetization 
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DOS vector at energy E reads 

™t, s (e) = £ £ g7T" ( *" 4 ' s ® )2/g2 / ^*§(^fe(^ (19) 

* ^ atomic volume 

with <r being the Pauli spin operator vector. Unit vectors determining the local spin quan- 
tization axis of the tip apex, er(-E'), and the ath surface atom, e^g(E), at a given energy can 
be calculated as 

eT,s{E) -=- - ^, (21)) 

m T ,s{E) x (£) 2 + m y (E) 2 + m z TS (E) 2 



thus, the atom-projected magnetization DOS vector can be rewritten as 



m T , s {E) = m T , s (E)e TtS (E) = yJm* TS (E) 2 + m y TS (Ef + m z TS (Efe T ^E). (21) 

Using this expression and the scalar product of the local spin quantization axes, 
eT(E)eg(E) = cosip a (E), the following holds, 

m T (E)m^(E) = m T (E)e T (E)m a s (E)e a s (E) = m T (E)m a s (E)cosy a (E). (22) 

The LDOS can also be written in terms of energy dependent spin polarizations. The 
spin polarization is defined as 

„ fJp . m TjS (E) nl s (E) - n l TS (E) 
PtAE) ~ ^AEj ~ 4 tS (E) + 4,sW' 
assuming collinear electronic structure. From noncollinear electronic structure the energy 
dependent spin polarization vectors are obtained by using Eq.(l2"T]) as 



m TS (E) J m T,s( E ) 2 + m T,s(E) 2 + "4,s(£) 2 

P T , S (E) = -^y- = ± — ^ e T , s (E) = P T AE)er,s{E). 

n T,s(ti) n TtS (E) 

(24) 
The relation between spin polarization vectors and scalars is similar to Eq. (l22[) . 

P T (E)P a s (E) = P T (E)e T (E)P^E)e a s (E) = P T {E)Pg(E)cos<p a (E), (25) 

with the same energy dependent unit vectors, which define the local spin quantization axes, 
see Eq.f l20p . Thus, the LDOS at the tip apex position and at energy E can alternatively be 
written using the above defined spin polarizations as 

LDOS{x,y,z,E,V) = /\E^e- 2K{Ey)da{x ^ z) n T {E)n a s {E)[l + P T {E)Pg{E)coscp a {E)}. 

a 

(26) 



Using Eq.(ll) of Ref. [14[ and our LDOS expression, a virtual differential conductance 
at the tip apex position and at energy E can be defined as 

^(x,y,z,E,V) = jY, e ~ ME ' V)da{X,V ' Z)nT W AEn ^ 

a 

(27) 
This means that by multiplying the LDOS with AE results in a dimensionless quantity, 
which is multiplied by the conductance quantum e 2 /h in order to arrive at our dl/dU 
expression. Note that ut(E)AE electron states from tip and n^(E)AE states from each 
surface atom contribute to the differential current at energy E, and in our model, dl/dU 
is proportional to the LDOS, which contains both surface and tip electronic information. 
If the two subsystems are calculated separately, it is possible to combine different levels 
of electronic structure for tip and surface, see also Table [H or include simplified model 
tip electronic structures into our approach. For example, assuming an electronically flat 
maximally spin-polarized (Pt{E) = 1) ideal magnetic tip with e.g. ht{E)AE = 1, the 
differential conductance reads 

^j(x,y,z,E) = j Y,^ ME)da{x ' y ' z) n a s {E)AE[l + P^{E)cos ( f a {E)]. (28) 



Here, Eq. ljTTj) has been assumed for the vacuum decay, and there is no ^-dependence. 

By measuring the energy with respect to the sample Fermi level as E = E F + eU, the 
energy dependence can be transformed to bias dependence U, as 

^-(x, y,z,U) = jJ2 e~ 2K ^ +eU ^ x ^ (29) 

a 

x n a s (E s F + eU)AE[l + P« {E S F + eU)cos^ a (E s F + eU)\. 
Similarly, the more general differential conductance, Eq.f l27|) . can be recast as 

^L(x,y,z,U,V) = e l^ E ) 2 Y J Z- 2K(ESF+eUy)da{x ' y ' z) (30) 

a 

x n T {El + eU- eV)n a s {E s F + eU)[l + P T (E^ + eU - eV)P%{E s F + eU)cos<p a (E F + eU)], 

where we used the fact that the tip Fermi level is shifted by eV with respect to the sample 
Fermi level, i.e. E F = E F + eV, and therefore E = E F + eU — eV. 

Virtual differential conductances have to be determined at E t points in a fine energy 
grid with AE resolution within an energy window [Ei(V, T), E 2 (V, T)} corresponding to the 
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applied bias voltage (V) and temperature (T). A value of 10~ 3 eV has been used for AE in 
our calculations, while we tested a finer grid (AE = 10 -4 eV) as well, with no improvement 
of our results. Finally, the tunneling current can be determined by the following energy 
integral, 

E 2 (V,T) 

I(x,y,z,V,T) = [ ^M, V) = AE y, ^ z ^ = 

J c du e *—' du 

Ei(V,T) E 1 <E i <E 2 

6 -(AEf J2 J2 e ~ MEi ' V)da(X ' V,Z) M E >s( E ^ 1 + WW)a^ a (^)], 

E 1 <E i <E 2 a 

where the energy window is defined as 

E X (V, T) = min (E%, E S F + eV) -In (3+ v 7 ^) k B T, (32) 

E 2 (V, T) = max (E F , E S F + eV) + In (3 + v^) A; B T. (33) 

Here, Ep is the Fermi energy of the sample surface, and V is the applied bias voltage. 
Broadening of electron states at finite temperatures is considered according to Eqs.( TT7l) . 
( TT8l) . ( fl9l) . and the temperature dependent terms in the integral limits are the full width 
at half maximum of the energy-derivative of the Fermi distribution function divided by 2, 
and ks is the Boltzmann constant. Another, more precise way to include thermal effects 



in calculating the tunneling current is given in the Appendix of Ref. 3l| based on the 

Sommerfeld expansion, which can also be incorporated into our approach. Lattice vibrations 

at nonzero temperatures are not taken into account. 

From the obtained three-dimensional tunneling current maps, data can be extracted 

which are directly comparable to experiments. In particular, current values can be shown 

in arbitrary z = Zq = const, planes or constant-value surface contours can be defined. The 

first option corresponds to the constant height mode, I(x,y,Z c = const. ,V,T), while the 

second to the constant current mode of SP-STM, I(x,y,z,V,T) = Ic = const. From the 

latter, a two-dimensional tip position map, called the height profile, z(x,y,V,T,Ic) can be 

extracted using logarithmic interpolation between grid points z\ < 22, if I(x,y,z\^ V,T) > 

Ic > I(x, y, z 2 , V, T), in the following way, 

/ w rp T x * ln(I c ) -ln{I{x,y,zi,V,T)) 

**.«M.'c) = * l + *z imx ^ VtT)) _ imXty ^ v ^ y (34) 

where Az = z 2 — Z\ = z i+ i — Zi = 0.0529177A (0.1 a.u.) has been used in all calculations. 
Alternatively, if Ic has such a value which is not contained in the considered finite box above 
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the surface, i.e. if Ic < I(z max ), then 

z(x, y, V, T,U- z max I Az ln{Ic) - HI(X ' V ' Zmax ' Yl T)) (35) 

ln(I(x, y, z max , V, T)) - ln(I(x, y, z max _ x , V, T)) 

Surface corrugation can be determined from this z(x,y) map. Note that the total current 
contains both to pog raphic and magnetic contributions, therefore z(x, y) will be the simulated 



SP-STM image 15] . in our model we consider its bias and temperature dependence as well. 
. n periodic magnetic systems the magnetic unit cell can be identified in the simulated image 
141 ] . Moreover, two different types of magnetic contrast can be defined. The first one is 
the apparent height difference of a particular atom at (xi,yi) lateral position imaged with 
a magnetic tip with parallel (P) and antiparallel (AP) relative magnetic orientation, on the 
same constant current contour, 

Azi(V, T, Ic) = z(x h Vu V, T, I p = I c ) - z(x h Vi , V, T, I AP = I c ). (36) 

The other magnetic contrast is the apparent height difference of two magnetic atoms at 
different lateral positions (x;, yi) and (xj, yj) imaged with a fixed tip magnetization direction, 

Azij(V, T, I c ) = z{xj,yj, V, T, I = I c ) - z(x h y h V,T,I = I c ). (37) 

This means apparent height differences of different surface atoms on the same SP-STM 



image, similarly as considered e.g. for oppositely magnetized islands in Ref. |32|. Note that 
both magnetic contrasts depend on the bias voltage, temperature and the constant current 
value. This latter means, in effect, dependence on the tip-sample distance. 

The obtained tunneling current can also be decomposed into a non-spin-polarized 
(TOPO) and a spin-polarized (MAGN) part, 

Itotal(x, y, z, V, T) = I T opo(x, y, z, V, T) + I MAGN (x, y, z, V, T), (38) 

with 

I TO po(x,y,z,V,T) = ^{AEf £ ^ e - 2 ^'^(^)n T (^)n°(E 4 ), (39) 

E 1 <E i <E 2 a 



T M AGN(x,y,z,V,T) = -(AE) 3 ]T ^e" 2K ^^^( a: *^m T ( J E; j )m°( J E; j ) C os^(^$40) 



E\<Ei<E2 a 

Itopo an d Imagn can be analyzed separately using the same way as described for the total 
current and they can be related to SP-STM experiments using the differential magnetic 
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mode J9|. From the non-magnetic height profile, z(x,y,V,T, Itopo = const.), the surface 
topography can be calculated, and in periodic systems the chemical unit cell is revealed in 
the simulated image. 

It has to be noted that the presented method can also be applied to study nonmag- 
netic systems, where all magnetic contributions are equal to zero and the corresponding 
topographic STM images can be simulated. 



Finally, it is important to note that following Ref. 3l|, the real physical differential 
conductance measured in experiments can be obtained as the derivative of the tunneling 
current, Eq. (|3ip . with respect to the bias voltage. It can be related to our virtual differential 
conductance dl/dU defined in Eg. (130)) in the following way, 

£(,, „, z, V, T = 0) = £(,, y, z, V, V) + [ CU A§(,, y, ,, V, V)\ v ^ . (41) 

Here, T = K temperature is considered for the reason of simplicity in the integral limits. 
In the case of assuming an ideal magnetic tip, i.e. if ut{E) = const., WLt(E) = const., 
and k(E) has no ^-dependence as defined in Eq. ljTTj) . then dl/dU has no ^-dependence 
as in Eq. (l29[) . and consequently, the integral term is zero. In that highly idealized setup, 
dI/dV(V) = dI/dU(V), such that Eq. (l29j) is the real differential conductance with U = 
V = V . Moreover, dl/dU can also be written as a sum of TOPO and MAGN parts, 

T ^j AL (x, y, z, U) = T °*° (x, y, z, U) + M J^ N (x, y, z, U), (42) 

with 

^|P(s, V, z, U) = jAE J2 e- 2 « E F +eU ^^n a s (E s F + cU), (43) 

a 

^f^(x, y, z, U) = 6 jAE J2 e-^ E ^ u ^>^m» s (E s F + eU)cos<p a {E% + eU)}, (44) 

a 

and they can be analyzed separately. We return to the simulation of spin-polarized scanning 
tunneling spectroscopy (SP-STS) based on the atom superposition method in the future. 

III. RESULTS AND DISCUSSION 

In order to demonstrate the capabilities of our model for simulating SP-STM on complex 
magnetic surfaces, we consider a sample surface with noncollinear magnetic order. One 
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ML Cr on Ag(lll) is a prototype of frustrated hexagonal antiferromagnets [151 ] . Due to 
the geometrical frustration of the antiferromagnetic exchange interactions between Cr spin 
moments, its magnetic ground state has been determined to be a noncollinear 120° Neel state 
[141 ] . We consider two possible Neel states with opposite chiralities, which are energetically 
equivalent only in the absence of spin-orbit coupling. 

We performed geometry relaxation and electronic structure calculations based on Density 
Functional Theory (DFT) within the Generalized Gradient Approximation (GGA) imple- 



mented in the Vienna Ab-initio Simulation Package (VASP) [34j-|36[. A plane wave basis set 
for electronic wavefunction expansion together with the projector augmented wave (PAW) 
method [37]] has been applied, while the exchange-correlation functional is parametrized ac- 



cording to Perdew and Wang (PW91) 38 



structure we used the VASP code as well 



39 



? or calculating the fully noncollinear electronic 
, |40j , with spin-orbit coupling considered. This 
allows us to determine the Neel state with the energetically favored chirality. 

We model the Cr/Ag(lll) system by a slab of a five-layer Ag substrate and one-one 
monolayer Cr films on each side, where the surface Cr layers and the first Ag layers under- 
neath have been fully relaxed. After relaxation the Cr-Ag interlayer distance is reduced by 
9.5%, while the underneath Ag-Ag increased by 0.5% compared to bulk Ag. A separating 
vacuum region of 14.6 A width in the surface normal (z) direction has been set up between 
neighboring supercell slabs. The average electron workfunction above the Cr is calculated to 



be 4> s = 4.47 eV using Eq.(E[3]). We used an 11 x 11 x 1 Monkhorst-Pack (MP) [4J| k-point 
grid for calculating the projected electron DOS onto the surface Cr atoms in our (v^3 x \/3) 
magnetic surface unit cell. 

Performing fully noncollinear electronic structure calculations we obtained convergence 
to two different magnetic Neel states. The magnetic surface unit cell with the converged 
magnetic moment directions are shown in the left part of Figure [TJ Each of the two Neel 



states can be characterized by a chirality vector, defined as 12| 



K = ^ {e l s x e| + e| x e| + e| x 4) . (45) 

Here eg denotes the local spin quantization unit vector of the ath Cr atom. It is defined 
from the local magnetic moment, M s = J_ F dEfng(E), similarly as in Eg. (1201 . i.e. e| = 
M S /\M s \. The magnitude of the magnetic moments of the Cr surface atoms are 3.73 yU^, 
with a very small out-of-plane component, which is neglected when defining the chirality 
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vectors. Thus, in the first row of Figure [TJ e\ = (1/2, y/%/2, 0), e| = (1/2, -y/%/2, 0), and 
e| = (—1, 0, 0). This corresponds to the chirality vector K = (0, 0, —1) or simply K z = —1. 
Similarly, in the second row of Figure dj e^ = (1/2,-^3/2,0), e| = (1/2, V3/2, 0), and 
e| = (—1, 0, 0) correspond to if z = +1. Comparing total energies of the two states we find 
that K z = — 1 is energetically favored by 1.1 meV compared to K z = +1. This finding 
is consistent with the magnetic ground state found for a Cr trimer island on the Au(lll) 



substrate in Ref. (12J, where it was also shown that the Dzyaloshinskii-Moriya interaction 
is responsible for determining the ground state magnetic chirality. Performing a collinear 
calculation with spin-orbit coupling considered, we obtain a ferromagnetic (FM) state with 
in-plane Cr atomic magnetic moments of 3.76 /i#. It turns out that this FM state is 1.04 eV 
higher in energy than the K z = — 1 Neel state. The energy difference of 346 meV/ (magnetic 
atom) in favor of the Neel state is in good agreement with results of Ref. [14J. The out- 
of-plane FM state is 1 meV higher in energy than the in-plane FM state with the same 
magnitude of magnetic moments. 

Simulation of SP-STM images can be performed using Eqs. (l27|) and ( 13~T|) in two ways: 
(1) According to Heinze 



15| . having chemically equivalent surface atoms the spin structure 
plays a more dominant role compared to the detailed electronic structure in determining 
the main features of an SP-STM image. Following this, we can take the collinear electronic 
structure (COLL) obtained from the in-plane ferromagnetic calculation, and set the spin 
structure to the corresponding Neel state. 

(2) As a more precise way, we can take the noncollinear electronic structure (NONCOLL), 
and there is no need to prescribe the spin structure as it is naturally included in the elec- 
tronic structure data. 

The first approach is computationally cheaper, and can be applied to simulate larger scale 
images 



15[. On the other hand, calculation of the noncollinear electronic structure is com- 
putationally more demanding but more realistic. 

In our Cr/Ag(lll) system we calculated the tunneling current in a box above the magnetic 
unit cell containing 153000 (34x30x150) grid points with a 0.15^4 lateral and 0.0529177A 
horizontal resolution. Figure [TJ shows simulated constant current SP-STM images for the 
two Neel states at zero bias voltage, assuming an ideal electronically flat maximally spin- 



polarized tip based on Eq. f[28|) with various magnetization directions following the first 
method. These are in qualitatively good agreement with previous simulations 14J, Il5| . 
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Using a nonmagnetic tip, all surface Cr atoms appear to be of equal height (one height 
level), i.e. the surface topography is seen. As the spin polarization of the Cr atoms at 
the Fermi energy is positive (see Figure |2]), and the tip spin polarization is set to +1, the 
Cr atom with parallel/antiparallel magnetization direction relative to the tip appears to 
be higher/lower than the other two Cr atoms, which have the same apparent height due 
to symmetry (two height levels). Comparing the images, it is clear that a contrast reversal 
occurs when turning the tip magnetization to opposite direction. This magnetic contrast can 
be quantified according to Eq.( )36|) . By setting the tip magnetization direction perpendicular 
to a Cr magnetic moment, a structure with three height levels occurs. This means that all Cr 
atoms in the magnetic unit cell have different apparent heights. This is due to the variation 
of the angles between the local Cr magnetic moments and the tip magnetization, e.g. for 
K z = —1, tpx = 90°, <f 2 — 30°, and <p3 = 150°. Determining the chirality of the magnetic 
structure from experimental SP-STM images is only possible in such a scenario if the tip 
magnetization direction is not parallel with the magnetic moment of any of the surface atoms. 
In our example of the Cr/Ag(lll) system the three apparent height levels follow a different 
order in the magnetic unit cell corresponding to the different chiralities. The decreasing 
levels of Cr apparent heights are indicated by circular arrows in the last column of Figure [TJ 
Apparent height differences of individual atoms on the same image define another kind of 
magnetic contrast, see Eq. fl37|) . Generally, the determining factor for the apparent height of 
magnetic atoms in zero bias (^=0 V) measurements is the effective spin polarization (ESP) 
at the common Fermi level, PT(Ep)Pg(Ep)cosip a (Ep), similarly as it was identified as the 



governing factor for the height of differential tunneling spectra at particular energies 24f . A 
positive ESP results in higher tunneling current at a fixed distance above a magnetic surface 
atom, while the opposite holds for negative ESP. Considering a constant current contour, 
thus, results in a higher apparent height for the atom with positive, while a lower height 
with negative ESP, compared to the topographic heights. 

Let us analyze the consequences of the choice of the collinear (COLL) or noncollinear 
(NONCOLL) electronic structure for the SP-STM images in more detail. Taking the non- 
collinear electronic structure we obtained qualitatively similar images at zero bias as shown 



in Figured], thus, Heinze is right [15| with the quality of the SP-STM images calculated at the 
sample Fermi energy using either COLL or NONCOLL electronic structure. The different 
spin polarization value of the Cr atoms at the Fermi level, however, results in different mag- 
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netic contrasts. According to Figure[2j the spin polarization of the Cr atoms is 0.20 and 0.51 
considering the COLL and NONCOLL electronic structure, respectively. Consequently, we 
expect that the magnetic contrast is higher in the NONCOLL SP-STM image. Indeed, e.g. 
taking a constant current contour of 10 -4 nA at parallel tip magnetization direction to the 
surface Cr atom labeled by "1" (Crl), we find Azg OLL (V = 0V, T = 4.2K, I c = 10" 4 nA) = 
0.07 A and Az^ 2 ONCOLL (V = 0V,T = 4.2K, I c = 10" 4 nA) = 0.21 A magnetic contrasts 
for COLL and NONCOLL images, respectively. Moreover, the 10 -4 nA contour is closer 
to the sample surface in the NONCOLL case. It is worth to compare magnetic contrasts 
of COLL and NONCOLL images on constant current contours having the same apparent 
height for Crl. For example, the Crl apparent height of 3.35 A is obtained at 10~ 4 nA in the 
COLL and 5 x 10 -5 nA in the NONCOLL image. The magnetic contrast in the new contour 
of the NONCOLL image is Az? 2 ONCOLL (V = 0V,T = 4.2K, J c = 5 x 10~ 5 nA) = 0.18 A. 
Thus, we find that the magnetic contrast ratio of NONCOLL and COLL images at the same 
Crl apparent height of 3.35 A, 0.18A/0.07A equals to the spin polarization ratio of Crl 
NONCOLL and COLL electronic structures at the Fermi level, i.e. 0.51/0.20. 

In the following we consider the magnetic Neel state with K z = — 1 chirality since it 
has been identified as the ground state. Figure |2] compares the energy dependent spin 
polarization vectors of Crl in Figure [IJ calculated from NONCOLL and COLL electronic 
structures. The spin polarization vector is defined as P S (E) = Pg(E)eg(E), where Ps{E) 
is calculated using Eq.( j23l and Eq. (j24l) in the COLL and NONCOLL case, respectively. 
Taking the collinear electronic structure, the local spin quantization axis of Crl is set to 
the local magnetic moment direction with neglecting the small out-of-plane component, 
e^ = (1/2, a/3/2, 0), and it is independent of energy. Reversal of the spin polarization 
vector occurs at Ps{E) values of opposite sign. Note that three sign changes occur in 
the [0.0 eV,0.3 eV] energy interval with respect to the Fermi level, using an 11 x 11 x 1 
MP k-point grid. We tested a denser 15 x 15 x 3 MP k-point grid as well, resulting in a 
qualitatively similar spin polarization. For computational and comparison reasons, we chose 
the llxllxl MP k-point set for calculating the NONCOLL electronic structure. While, in 
the noncollinear case P${E) is always positive due to Eg. (1241) . one spin polarization vector 
reversal is observed at 0.54 eV above the Fermi level. The indication for this reversal is 
the sign change of e s x (E) and eg(E), i.e. going away from the Fermi energy, the local spin 
quantization axis changes from e^ ~ (1/2,-^/3/2,0) to e^ ~ (— 1/2, — \/3/2, 0) at 0.54 eV. 
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Here, however, the e s z (E) components are not exactly zero, but they are in the order of 
10 -6 to 10~ 2 in the whole energy range with the exception of e<f (0.54eV) = —0.21. This 
latter value indicates that the rotation direction of the spin polarization vector at 0.54 eV is 
through negative e s z components. Since e s (E) is a unit vector at all energies, the presence 
of small e s z (E) components also means that the other vector components are |e;!f (_E)| » 1/2 
and \eg(E)\ ~ v3/2. By comparing the Pg(E) spin polarization function of COLL and 
NONCOLL electronic structures, we can state qualitative agreement. 

Let us compare simulated single point differential conductance spectra based on NON- 
COLL and COLL electronic structures. Figure [3] shows such simulated spectra z = 3.5 A 
above the Crl atom in FigureQ]with assumed parallel (P) and antiparallel (AP) tip magneti- 
zation direction using an ideal magnetic tip. We showed at the end of section [Til that for the 
considered ideal magnetic tip dl/dU (Eg. ( 129]) ) is the real differential conductance. Accord- 
ing to Eq. ()4"2]) . the topographic and magnetic contributions can be calculated separately. 
Determining dlTOPo/dU(z,U) (red dashed line with symbol "X") and dl[ iAGN /dU(z,U) 
(blue dashed line with symbol "+") is sufficient to draw dl p /dU(z, U) (black solid line) and 
dI AP /dU(z, U) (brown (gray) solid line) since 

dI p /dU( Zl U) = dI TO Po/dU(z, U) + dI^ AGN /dU(z, U) 
dI AP /dU(z, U) = dI TOPO /dU(z, U) - dI p AGN /dU(z, U). (46) 

Here, we took into account that the magnetic contribution for the AP tip magnetization di- 
rection dI A f AGN /dU equals to —dI^j AGN /dU, since cosip changes sign. The COLL and NON- 
COLL spectra have slightly different peak positions due to the details of the electronic struc- 
ture. We find that dI p /dU > dI TO Po/dU > dI AP /dU below U = 0.54 V, while dI p /dU < 
dlroPo/dU < dI AP /dU above U = 0.54 V, calculated by using NONCOLL electronic struc- 
ture. The relation of these quantities is determined be the sign of the magnetic contribution 
at the given bias, i.e. dI^ AGN /dU(U < 0.54V) > 0, and dI^ IAGN /dU(U > 0.54V) < 0. On 
the other hand, for the COLL case, there are three sign changes of dI p IAGN /dU at 85 mV, 
160 mV, and 300 mV, resulting in dI p /dU > dI TO Po/dU > dI AP /dU below U = 85 mV 
and dl p /dll < dlxopo I dU < dI AP /dU above U = 300 mV. The magnetic contribution is 
small between 85 mV and 300 mV and the difference between spectra is less than 0.02 nA/V 
in this bias range. 

In the following we use the NONCOLL electronic structure for the Cr/Ag(lll) sample 
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surface. By including energy dependent electronic structure of sample and tip into our 
model, we can study the bias dependent magnetic contrast and its tip dependence as well. 
Figure H] shows simulated SP-STM images for various tip magnetization directions at -1 V, 
V and +1 V bias voltages assuming an ideal magnetic tip. We find qualitatively similar 
images for -IV and V for the corresponding tip magnetization direction. This means 
that the Azi2 magnetic contrast between Crl and Cr2 has the same sign at -1 V and V. 
However, Azi2 at the same Crl apparent height increases at -1 V compared to V for all 
tip magnetization directions. This can be explained by the integrated dI^ IAGN contribution, 
which does not change sign in this bias range, see blue dashed line with symbol "+" in the 
left part of Figure On the other hand, the results show that the magnetic contrast is 
reversed at +1 V compared to the other two studied bias voltages. This contrast reversal 
is observed for all tip magnetization directions. It is interesting to find that on the image 
with three height levels the apparent heights change order in such a way that the image at 
+ 1 V looks like that the Neel state would have an opposite chirality compared to V or 
-IV, see the indicated circular arrows in the last row of Figure HI This finding highlights 
the importance of the applied bias voltage and suggests that one has to be careful when 
interpreting the magnetic structure from experimentally observed SP-STM images. Based 
on our theoretical study we can also conclude that the magnetic contrast reversal occurs 
between V and +1 V bias voltages. This contrast reversal is solely due to the sample 
electronic structure since the ideal magnetic tip is electronically featureless. 

Dependence of the magnetic contrast on the tip electronic structure can be studied by 
considering different tip models. As an example we chose a ferromagnetic Ni tip. Such tips 



are routinely used in SP-STM and STS experiments [42|, |43j . The Ni tip has been modeled 
by a seven-layer Ni film slab with (110) orientation, having one-one Ni apex atoms on both 
surfaces, i.e. with a double vacuum boundary. Here, the apex atom and the topmost surface 
layers have been relaxed on both sides. The interaction between apex atoms in neighboring 
supercells is minimized by choosing a 3 x 3 surface cell, and a 15.4 A wide separating vacuum 
region in z direction. Moreover, an llxllxl MP k-point grid has been chosen for obtaining 
the projected DOS onto the apex atom. The electronic structure of the apex is given in the 



top part of Figure 1 of Ref. [24|. We obtain a spin polarization of Pt = —0.91 at the Fermi 
level, E%, and \P T {E)\ > 0.8 between E T F - 0.3eV and E T F + 0.3eV. Employing Eq.flU}, the 
local electron workfunction above the tip apex is <pT = 4.52 eV, and Eq.( iT2l has been used 

19 



to determine the vacuum decay. 

Figure [5] shows simulated SP-STM images for various tip magnetization directions at 
-1 V, V and +1 V bias voltages including the electronic structure of the Ni tip into 
our model. By comparing images to those shown in Figure H] obtained by using an ideal 
magnetic tip, we find that the magnetic contrast is the opposite for each picture. This 
is due to the negative spin polarization (-0.91) of the Ni tip apex at its Fermi level 24]. 
Note that the spin polarization of the ideal magnetic tip was assumed to be +1 in the 
whole energy range. Similarly as in Figure HJ we find qualitatively similar images for -1 V 
and V for the corresponding tip magnetization directions with higher magnetic contrast 
at -1 V compared to V. Again, the magnetic contrast is reversed at +1 V compared to 
the other two studied bias voltages. This effect is highlighted in the last row of Figure 
showing the decreasing levels of Cr apparent heights by circular arrows, thus, indicating 
a bias dependent apparent magnetic chirality. The results suggest that different tips can 
completely reverse the magnetic contrast. This effect has to be taken into account when 
determining the magnetic structure from experimentally observed SP-STM images. 

IV. CONCLUSIONS 

We extended the atom-superposition-based method of Heinze [15j for simulating spin- 
polarized scanning tunneling microscopy by including the tip electronic structure, bias volt- 
age, and the capability of incorporating the fully noncollinear electronic structure. Taking 
the tip electronic structure into account, the effect of a richer variety of electronic structure 
properties can be investigated on the tunneling transport within the indicated approxima- 
tions (atom superposition, spherical vacuum decay). The method is computationally cheap 
and it can be applied based on results of any ab initio electronic structure code. Taking 
the prototype frustrated hexagonal antiferromagnetic system, Cr monolayer on Ag(lll) in 
a noncollinear magnetic 120° Neel state, we determined its ground state magnetic chirality 
and simulated SP-STM images at different bias voltages to illustrate the applicability of our 
method. We related the magnetic contrast of the zero bias images to the effective spin po- 
larization at the sample Fermi level. Moreover, we illustrated the importance of the energy 
dependent local spin quantization axes by comparing collinear and noncollinear electronic 
structure of a particular surface Cr atom and its effect on single point tunneling spectra. 
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Finally, we showed evidence that the magnetic contrast is sensitive to the tip electronic 
structure, and this contrast can be reversed depending on the bias voltage. 
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TABLE I: Combinations of taking into account energy dependence (Y) or independence (N) of the 
local spin quantization axes of surface atoms and the tip apex, and consequence for the energy 
dependence of the angle <p a between the surface local and the tip quantization axes. The combi- 
nations framed by black solid lines are considered in the present work, while the one with the gray 



dashed line corresponds to the studied system in Ref. 
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FIG. 1: (Color online) Surface geometry of 1 ML Cr on Ag(lll) and simulated SP-STM images at 
V bias voltage depending on the tip magnetization direction (M T j P ) assuming an ideal electron- 
ically flat maximally spin-polarized tip. The Cr and Ag atoms are denoted by spheres colored by 
green (medium gray) and purple (dark gray), respectively, while the magnetic moments of individ- 
ual Cr atoms are indicated by (red) arrows in the left part of the figure. The Cr atoms are explicitly 
labeled corresponding to the calculated chirality vector in Eq. (|45p . In addition, the (v3 x y/3) 
magnetic unit cell is drawn by yellow (light gray) color. In the two rows noncollinear Neel states 
with opposite chiralities and corresponding SP-STM images are shown. In the last column, the 
decreasing levels of Cr apparent heights are indicated by circular arrows. 
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FIG. 2: (Color online) Comparison of energy dependent spin polarization vectors, P S (E) = 
P^(E)(e 1 s x (E),e 1 s y (E),e 1 s z {E)), of the surface Cr atom labeled by "1" in Figure [Q in the Neel 
state with K z = — 1 chirality, calculated from noncollinear (NONCOLL) and collinear (COLL) 
electronic structure. In the collinear case the local spin quantization axis is e^ = (1/2, \/3/2, 0) 
in the basis of (e x ,e y ,e z ), and is independent of energy. Reversal of the spin polarization vector 
occurs at P${E) values of opposite sign. In the noncollinear case Ps{E) is always positive due 
to Eq. (|24p and the spin polarization vector reversal is observed as the sign change of e l g{E) and 
eg(E) at 0.54 V. Here, the rotation direction of the spin polarization vector is through negative 
eg z (E) components. 
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FIG. 3: (Color online) Simulated single point differential tunneling spectra dl /dU and dl jdU 
3.5 A above the surface Cr atom labeled by "1" in Figure Q] in the Neel state with K z = — 1 
chirality, assuming parallel (P) and antiparallel (AP) tip magnetization direction with respect 
to that of Crl, applying an ideal electronically flat maximally spin-polarized tip according to 
Eq. (|29p . Left and right parts correspond to spectra obtained from noncollinear (NONCOLL) and 
collinear (COLL) electronic structures of the sample, respectively. Topographic (dlroPo/dU) 
and magnetic (dlf IAGN /dU = —dI^ AGN /dU) contributions are given according to Eq. (jl3|) and 
Eq. (j44|) . respectively. 
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FIG. 4: (Color online) Simulated SP-STM images depending on the bias voltage and the tip 
magnetization direction assuming an ideal electronically flat maximally spin-polarized tip. The 
magnetic contrast is reversed between 0.0 V and 1.0 V. In the last row, the decreasing levels of Cr 
apparent heights are indicated by circular arrows. The surface geometry of 1 ML Cr on Ag(lll), 
its magnetic structure with K z = — 1 chirality, and the considered tip magnetization directions are 
explicitly shown, similarly as in Figured) 
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FIG. 5: (Color online) Simulated SP-STM images depending on the bias voltage and the tip 
magnetization direction assuming a model Ni tip. The magnetic contrast is reversed compared 
to images obtained by using the ideal magnetic tip (compare to Figure ED, and there is a bias 
dependent contrast reversal between 0.0 V and 1.0 V. In the last row, the decreasing levels of Cr 
apparent heights are indicated by circular arrows. The surface geometry of 1 ML Cr on Ag(lll), 
its magnetic structure with K z = — 1 chirality, and the considered tip magnetization directions are 
explicitly shown, similarly as in Figure [TJ 
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